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It was previously suggested that an odd-frequency pair amplitude exists in the vicinity of 
boundaries in unconventional superconductors. We develop this idea and quest for a novel 
superconducting order parameter with an odd-frequency dependence. For this purpose, we focus 
on p-wave superconductors and extend the quasi-classical theory to include the odd-frequency 
dependence in the order parameter. Both of the frequency and spacial dependences of the 
order parameter are determined self-consistently. Under a finite electron-phonon interaction, it 
is found that an odd-frequency order parameter is stabilized near the boundary and coexists 
with the even-frequency one. By analyzing the induced odd-frequency pair amplitude in terms 
of the superconducting quasi-particle wavefunction, it is found that the mid-gap bound state 
generates the emergent odd-frequency order parameter. 

KEYWORDS: odd-frequency superconductivity, unconventional superconductivity, chiral p-wave supercon- 
ductivity, spin triplet s-wave, quasi-classical Green's function, Andreev equation, mid-gap 
state 



1. Introduction 

Odd-frcqucncy pairing state was suggested by Berczin- 
skii as a possible candidate to describe the superfluid 
•^He.^ Although it was not adapted to the superfluid '^He, 
the idea was developed and applied also to conduction 
electron systems. ^"^'^ In these works, the main subject of 
interest was to understand how the odd-frequency super- 
conductivity is realized in the bulk systems. 

From a different point of view, Bergeret et al. sug- 
gested that an odd-frequency pair amplitude is induced 
near the interface between superconductor and ferromag- 
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Since the magnetic field from the ferromagnet 
modifies the spin-part wavefunction of the Cooper pair, it 
leads to a mixing between the singlet and triplet parings 
and induces the odd-frequency pair amplitude. It was 
also pointed out by Tanaka et al. that the odd-frequency 
amplitude is present also in the vicinity of boundaries in 
unconventional superconductors. In this case, the 
even- and odd-parities of the Cooper pair are mixed by 
the broken translational symmetry. This parity mixing is 
also present around a vortex core and induces the odd- 
frequency amplitude there. ^^'^^ 

In the previous works, they focused only on the odd- 
frequency pair amplitude and did not pay attention to an 
order parameter. To understand this, it is important to 
distinguish between the pair amplitude and the order pa- 
rameter explicitly. The pair amplitude (F) is defined as 
an expectation value of the field operators of the Cooper 
pair and is understood as a superconducting correlation 
function, while the order parameter (A) is defined as a 
product of the pair amplitude and a coupling constant 
{V), i.e. A = VF{A). Since the pair amplitude is a func- 
tion of A, the order parameter should be determined in 
a self-consistent calculation. 

* E-mail address: spmmatu@ipc.shizuoka.ac.jp 



In the presence of the pair amplitude, a supercon- 
ducting order parameter is stabilized under a finite in- 
teraction between conduction electrons. In fact, we in- 
vestigated an electron-phonon mediated s-wave super- 
conductor under a finite magnetic field on the basis of 
the Eliashberg theory and found that an odd-frequency 
spin-triplet s-wave order parameter coexists with that 
of an even-frequency spin-singlet s-wave. ^^'^^ The mag- 
netic field leads to the singlet and triplet mixing and 
induces the odd-frequency pair amplitude. Since the 
electron-phonon interaction leads to an attractive inter- 
action also in the spin-triplet s-wave channel, the odd- 
frequency order parameter is stabilized. In other words, 
it can be understood as a consequence of the even- and 
odd-frequency mixing of the order parameter under the 
broken time-reversal symmetry.^^'^^ Physical quantities, 
such as transition temperature, density of states, re- 
sponse to external fields, show a different behavior owing 
to the odd-frequency order parameter. 

Similarly to the bulk case, we propose in this paper 
that the even- and odd-frequency order parameters co- 
exist in the vicinity of boundaries or interfaces of un- 
conventional superconductors. To demonstrate this, we 
focus on surface and domain wall in two-dimensional p- 
wave superconductors and study how the odd-frequency 
order parameter appears. To show this, we treat the p- 
wave order parameter on the basis of a weak coupling 
theory assuming that the p-wave order parameter has 
no frequency dependence. On the other hand, we have 
to take the frequency dependence into account in the 
odd-frequency order parameter. It is known that the sur- 
face induces an odd-frequency spin-triplet s-wave pair 
amplitude in the p-wavc superconductors.^'^ As the at- 
tractive interaction of the spin-triplet s-wave, we con- 
sider an electron-phonon interaction as in our previous 
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study under an external magnetic field. We deter- 
mine both the frequency and spatial dependences of the 
odd-frequency order parameter self-consistently on the 
basis of our previous quasi-classical theory for the sur- 
face and domain wall in the p-wave superconductors.^^ 

To understand the emergence of the odd-frequency or- 
der parameter, we express the induced odd-frequency 
pair amplitude with the superconducting quasi-particle 
wavefunction. It enables us to see what type of energy 
eigenstate generates the odd-frequency pair amplitude. 

This paper is organized as follows. In Sect. 2, the quasi- 
classical formulation for the odd-frequency order pa- 
rameter is presented. The solution of the quasi-classical 
Green's function is shown in Sect. 3. We discuss the ori- 
gin of the emergent odd-frequency order parameter in 
Sect. 4. The last section gives summary and discussions. 

2. Formulation 

2. 1 Odd-frequency pair amplitude 

Let us consider a two-dimensional p^j-wave supercon- 
ductor in a; > region and assume that a specular sur- 
face along the y-direction is located at a: = 0. Owing 
to the translational symmetry along the y-direction, the 
order parameter depends only on x. In this geometry, 
it is known that there is a mid-gap bound state in the 
superconducting energy gap near the surface. The quasi- 
classical Green's function shows this point clearly. For a 
d-vector parallel to the z-axis, the quasi-classical Green's 
function is separated into two identical 2x2 matrices. Al- 
though the self-consistent p^j-wavc order parameter has a 
spatial dependence, we assume a uniform one to see the 
induced odd-frequency pair amplitude in a simple way. 
For a constant order parameter, the spatial dependent 
quasi-classical Green's function is expressed as^^"'^^ 



E 

1=1,2,3 



9i{k 



F, 



(2.1) 



where pi represents the Pauli matrix of the i{= 1,2,3)- 
th component in the Nambu space, gi^kp, (jJm, x) is given 

by 

giikpx, kYy,UJrn,x) = gi{-kFx,kFy,UJ,n,x) 



|A.^(fci 



-2qx 



g2{kFx, kFy,UJm,x) = -.g2(-^Fa;, ^Fy, Wm,a;) 



-2qx\ 



A2(fcj 



gsikpx, fcFy, W„i, x) ^ g3{-kFx,kFy,UJ,n, x) 

, A^(fcF) 



^COm^+AUkF) [ 



-2qx 



,(2.2) 



with 



(2.3) 

Here, — 7rT(27Ti-f 1) (m: integer) is the fermionic Mat- 
subara frequency at temperature T . We take h = \ and 
fee = 1 throughout this paper, vfx represent the x com- 
ponent of the Fermi velocity, fcp = (kFx^ kpy) is the Fermi 
wave vector and its dependence in the pa;-wave order pa- 



rameter is expressed by Ax{k-p). It is characterized by the 
following symmetry: A^ikFx.kpy) = A^ikp^.-kFy) = 
—Axi—kpx, kpy)- We can see that Eq. (2.1) recovers the 
bulk solution when we take x — cx) as 



5(fcF,ii^m,bulk) 



y/cOrJ+Alikp) \-iA,{kF 



iA^(fcF) 



(2.4) 



The mid-gap {E = 0) bound state is described by the 
second term in 93 in Eq. (2.2). It decreases exponentially 
as e~^*^ for a: — > 00. Another characteristic point is that 
^2 = at .T = 0. This indicates that the Pa;-wave order 
parameter is destroyed by the scattering at the surface 
and vanishes at a; = 0. 

In addition to these, we can see that a finite gi compo- 
nent appears in the vicinity of the surface. According to 
the symmetry with respect to k^ —k^ and ky o — fcj,, 
the gi component has an s-wave symmetry. It indicates 
that the s-wave pair amplitude is induced near the sur- 
face. This s-wave pair amplitude is a spin-triplet one and 
has an odd symmetry with respect to Wm ^ —^m- Thus, 
the odd-frequency spin-triplet s-wave pair amplitude is 
induced near the surface. 

2.2 Odd- frequency order parameter 

As the typical unconventional superconducting state, 
we focus on two-dimensional p^,- and chiral Px + ipy-wave 
states, where the latter is suggested in Sr2Ru04 super- 
conductor. ^^"'^^ For simplicity, we assume a single cylin- 
drical Fermi surface. To determine the frequency and 
spatial dependences of the order parameter, we use the 
quasi-classical Green's function introduced by Schopohl 
et al. to study the vortex problem. We adapt this 
scheme to the boundary problem and apply it to both 
the cases of surface and the domain wall by extending 
our previous formulation^^ to include the odd-frequency 
order parameter. 

The momentum dependent p-wave state is given by 
the d-vector. In the bulk region, we assume a case of 
d(fc) = (0,0,Abuik(fc)) with Abuik(fe) = A^{k) + iAy{k), 
where A^(fc) and Ay{k) are order parameters for the 
Px- and Py-wave states, respectively. We take real values 
for Ax{k) and Ay{k). In the bulk region, the Matsubara 
Green's function is given by the following 4x4 matrix 
form: 



G(uJr,i,k) 



ekps - Abuik(fe)] ^ . 



(2.5) 



Here, Abuik(fc) is the matrix for the order parameter 
defined as 



Abuik(fc) = A^(fc)pi(Ti - Ay{k)p2(Ti, 



(2.6) 



where cr^ and pi {i = 1,2,3) are Pauli matrices for spin 
and Nambu spaces, respectively. For the p^: + ipy-wave, 
we apply weak coupling theory assuming that there is no 
frequency dependence. 

Let us discuss an additional order parameter appear- 
ing near surface. In this section, we consider a surface 
along the y-direction. In this case, the Pa;-wave order pa- 
rameter is suppressed near the surface owing to the odd 
parity with respect to the refiection at the surface. Since 
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the spin-triplet s-wavc pair amplitude is induced, we in- 
troduce the following the additional order parameter for 
the odd- frequency component: 



Aodd(fc,l^m) = As{k,ujm)ip2cri. 



(2.7) 



Here, As(A;,a;„i) is a real value. Owing to the s-wave 
symmetry, we assume that there is no k dependence in 
As(fc,aJm). For the odd-frequency order parameter, it is 
important to retain the frequency dependence. The s- 
wave order parameter satisfies 



As(fc,-w,„) = -As(A;,a;„i) 



(2.8) 



reflecting the odd-frequency nature. The matrix form of 
the order parameter is then given by 



A(fc,W„i) = Abulk(fe) + Aodd(fc,Wm) 

Ai(fc,w„i)(Ti 



A2(fc,a;„i)(Ti 







(2.9) 



There is the following relation between the order param- 
eters: 

A2(fc,cj™) = -At(-fc,cj™). (2.10) 

This conventional relation holds even when the odd- 
frequency order parameter is induced by the even- 
frequency one and exists as a minority component of the 
order parameter. ^^'^^ On the other hand, the relation in 
Eq. (2.10) does not hold anymore if the odd-frequency 
order parameter is the majority component. We have to 
use the other correct form in the latter case.^^'^^"^* Since 
the present problem corresponds to the former case, the 
order parameters are given by 

Ai(fc J ^m) — Ax{k) + iAy (fc) + As(fc, iOfn), 

A2(fe,w,„) = A:,(fc)-iAj,(fc)- A,(fc,cj„). (2.11) 
There are the following relations: 

Ai(-fc,-w„j) = -Ai(fc,i:j„), 

A2(-fc,-W„j) = -A2{k,UJrn), 

A2{k^,-ky,uj„i) = A2{kx,ky,uJrn)- (2-12) 

2.3 Quasi-classical theory with frequency dependent or- 
der parameter 
Since the odd-frequency order parameter appears only 
near the boundary, we have to treat a non-uniform situ- 
ation. For this problem, it is convenient to use the quasi- 
classical Green's function. For the order parameter given 
by Eq. (2.11), we obtain the following Eilenberger equa- 
tion: 



ivp ■ \/g{kF,ujm,r) 

iwm -Ai(fcF,w,„,r) 
A2{kF,ujrn,r) -iw™ 



,g{kF,LL>.m,r) 



(2.13) 

where g(k-p,uj„nr) is the quasi-classical Green's func- 
tion in a 2 X 2 matrix form.^^'^^ The kp dependence in 
Ai(feF,a;m,r) and A2(fcF,Wm,f") represents the orbital 



symmetry of the order parameter, while r in Eq. (2.13) 
represents the center of mass coordinate of the Cooper 
pair. The ujm dependence in the order parameter is owing 
to the odd-frequency component. 

The quasi-classical Green's function can be written as 



g{kF,u;rn,r) 



g{kF,uj,n,r) i/(fcF 
-i./(feF,Wm,r) -g(fcF,w™,r) 



(2.14) 



where the components (7(fcF, Wm, t), /(fcF,Wm,r), and 
f{k-p,ujrmr) satisfy the following equations: 

vp ■ \7g{kF,uj,n,r) = A2(fcF, t^m, r-)/(A;F, w,„, r) 
- Ai(fcF r-)/(fcF 

+ ]^VF ■ f{kF,U},n,r) 

= Ai{kF,ujm,r)g{kF,ujm,r), 



Wr„ - -vf ■ V ) f{kF,u;m,r) 



= A2{kF,uj,n,r)g{kF,ui,n,r). (2.15) 
They can be written as^^ 

1 - a{kF,LL>m,r)b{kF,LL>„i,r) 



g{kF,LOm,r) 
f{kF,ujm,r) 
J{kF,uj,n,r) = 



1 + a{kF,uJm,r)b{kF,uJni,r)' 

2a{kF,ujjn,r) 

1 + a(fcF, w™, r)b{kF,ujrn, r) 
2b{kF,uJm,r) 



(2.16) 



1 + a(feF,w„i,r)6(fcF,w,„,r) ' 
where a and b satisfy the following equations: 
vp ■ Va(fcF,w™,r) = Ai{kF,ujjn,r) 

- A2{kF,i^r,i,r)a'^{kF,u)m,r) - 2iu™a(fcF, w™, r), 
vp ■ Vfe(fcF, i^m, r) = -A2(A;f, Wm, r) 

-\- Ai(fcF, Wm, r)&^(fcF, Wm, r) -f 2cj™6(fcF, ujm, r). 

(2.17) 

We solve Eq. (2.17) along the quasi-classical trajectory 
as shown in Fig. 1(a), where the quasiparticle moves from 
A to B with the momentum kpi and from B to C with 
momentum A;f2 • The Fermi surface on the k^ — ky plane 
is cylindrical as shown in Fig. 1(b). Owing to the transla- 
tional symmetry along the y-direction, kpy is conserved. 
Then, we match the two solutions using the following 
boundary condition at point B:**^"^" 

g{kFl,UJ,n,B) = g{kF2,^rn,B), (2.18) 

which means for a and b as 

a{kFi,UJm,B) = a(feF2,'^m,B), 

6(fcFi,a;,„,B) = fe(fcF2,w™,B). (2.19) 

Since the system is translationally invariant along 
the y-direction, the quasi-classical Green's function or 
a{kF,uJm,r) and 6(fcF,Ct'm,r) depend only on x. Then, 
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Fermi surface 




(b) 



Fig. 1. (Color online) (a) Quasi-classical trajectory of a quasi- 
particle, in which the momentum of the incident and reflected 
quasiparticles along the surface is conserved, k-pi and fcF2 are 
the momentum with kp^<0 and kpj.>0, respectively. A specu- 
lar surface running along the y-direction is located at a: = 0. 
(b) Fermi surface and Fermi wavevector for the two-dimensional 
system. Sj. is measured from the A;a;-axis. 



Eq. (2.17) can be rewritten as 
ax 

~ A2{kF,uj„i,x)a'^{kp,uj,n,x) - 2u)„ia{kp,uj,n, x), 

VFx-^b{kF,UJm,x) = -A2(A;F,ii^m,a-) 

+ Ai(fcF,Wm,a;)6^(fcF,Wm,x) -|- 2uJmb{kF , UJm, x) . 

(2.20) 

where vfx is the x component of the Fermi velocity. The 
initial and boundary conditions for Eq. (2.20) are given 

by 



a{kFi,UJm, X = oo) 



Ai(fcFl,Wm,Oo) 



^^^2 + Ai(feFi,W™,CX))A2(fcFl,a;m,00) +UJm' 
h{kF2-,^m-,X = oo) 

A2(fcF2,'^rn,Oo) 



Ai(feF2,'^m,0o)A2(fcF2,W™,00) + U}„ 



2VxCosdkRe[f{ek,i 
2VySm0klTn [f{ek,L 



^,X) - f{TT - 9k,UJm,x)] 
i,x) + /(tT - 9k,UJm,x)] 



.VsV-{uji,uJm)'R.e [f{9k,i^m,x) + /(tt 



9fc,Wm,a;)] 

(2.23) 



where the following relations were used to derive the gap 
equation: 



f{9k + TT, -CJ, 
/( 



m,x) = -f{9k,UJ.m,x), 
f*{9k,LL>m,x). 



(2.24) 



In Eq. (2.23), Vx, Vy, and Vs are dimensionless coupling 
constants for the px-, Py- and s-waves, respectively. In 
the Px + ipy-wave pairing, the coupling constants are ex- 
pressed as^2 

-1 

E 17^1 ' (2-25) 

0<m<^ 



Vx^Vy = log 



T 



where Tc is the superconducting transition temperature 
for the p-wave. 

Next, we consider an effective interaction for the s- 
wave. Since the s-wave order parameter depends on the 
frequency, we need a frequency dependence in the in- 
teraction. As the simplest model, we consider Einstein 
phonons coupled to the conduction electrons. Since the 
electron-phonon interaction does not depend on spins of 
the conduction electrons, it is attractive for both the sin- 
glet and triplet channels. "'^^ The odd-frequency depen- 
dence of the effective interaction for the s-wave is given 
by"' 27 



(2.26) 



a(fcFi, Wm, a; = 0) = a(fcF2, ^m, a; = 0), 

b{kFi,UJm, X ^0) = 6(fcF2, I^m, X = 0). (2-21) 

We write the superconducting order parameters as 

Ai(fcF,a;,„, a;) = Ax{x) cos9k + iAy{x) sm9k + As(w,„, a;), 

A2(fcF,w„,a;) = A2;(a;)cos6'fc - iA^ (a;) sin 6*^ - As{uijn,x), 

(2.22) 

where A^ix), Ay{x), and As{u}m,x) are the supercon- 
ducting order parameters for the Px--, Py-, and s-waves, 
respectively. The Fermi wavevector dependences are ex- 
pressed by the angle of 9k as shown in Fig. 1(b). The 
order parameters are real numbers and are determined 
by 51 



where we is the frequency of the Einstein phonon. 
We note that V-{uJi,ujm) is antisymmetric with respect 
to LOi -uji and w™ -w™, i.e. y_(-w/,Wm) = 
V-{uJi,~ujm) = -V-{ui,Lo,n)- In Eq. (2.23), we intro- 
duced a cutoff energy Wc in the Matsubara frequency 
summation. In general cases, the cutoff energies are dif- 
ferent for the p- and s-waves, however, the cutoff energies 
are irrelevant to discuss the odd-frequency order param- 
eter. We assume the same cutoff energies for both the p- 
and s-wave order parameters here. 

Now we can determine the order parameters by solv- 
ing Eq. (2.23) with Eqs. (2.16), (2.20), and (2.21) self- 
consistently. 

3. Self-Consistent Solution 

3.1 p^-wave case 

First, we study a case of px-wave assuming Vy = 
[Ay{x) = 0]. Before going the non-uniform solution, let 
us summarize the bulk solution. The transition tempera- 
ture Tc and the solution of the gap equation [Eq. (2.23)] 
at T = for x — >■ oo are obtained as^^''^^ 

2(PuJc 1 
Ir = e , 



Axix) 
Ay{x) 

yAs{uJl,x) 



2T 



d9k 



Ax{oo) = AuJce 



2.Un = An. 



(3.1) 
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Fig. 3. (Color online) (a) Temperature dependence of the or- 
der parameters scaled by Aq. A^; is the pa;-wave component in 
the bulk region [x — > oo). As is the maximum magnitude of 
\As((^m,x)\ at x = 0, where the maximum value appears at 
f^m ~ 17Tc. (b) Coupling constant (Vs) dependence of the max- 
imum value of \As{uJm, ^)\ at x = for T = O.lTc. 



100 




Fig. 2. (Color online) (a) Spatial dependence of the p^-wave or- 
der parameter at T = O.lTc for Vy = (pure px-wave). The 
order parameters are normalized by the balk value at T = 0, 
Ao = Aa;(oo) ~ 2.14Tc. The x coordinate is measured in unit of 
the coherence length defined by § = iip/Ag. The parameter is 
chosen as ojc = lOOTc. (b) Spatial dependence of the s-wave order 
parameter. The Matsubara frequency is fixed as uJm = 17.3Tc at 
which \Ag{u}jn, x)\ has the maximum value. The set of param- 
eters are chosen as los = lOTc and Vs = lOVx- (c) Matsubara 
frequency dependence of the odd-frequency spin-triplet s-wave 
order parameter at various positions of x. (d) Contour plot of 

As{uJm,x). 



Here, 7 is the Euler's constant: 7=0.57721- •• . For the 
characteristic length, we define the following coherence 
length for the superconducting state: 

In Fig. 2, we show the self-consistently determined or- 
der parameters. The order parameters and the spatial 
coordinate arc scaled by Aq and ^, respectively. The set 



(3.2) 



of parameters are chosen as lJc = lOOTc, lue = lOTc, and 
Vs = lOVx- We can see in Fig. 2(a) that Aa:{x) decreases 
near the surface and disappears completely at x = 0.^^ 
In contrast to this component, the odd- frequency s-wave 
order parameter As(a;m, x) appears in the vicinity of the 
surface as shown in Fig. 2(b). Its magnitude has a maxi- 
mum value at a; = and decreases with x monotonically. 
In Fig. 2(c), we show the Matsubara frequency depen- 
dence of the s-wave order parameter at various positions 
of X. Since the s-wave order parameter must have an 
odd-frequency dependence owing to the fermion prop- 
erty. As {uJm , x) shows a linear Wm dependence for small 
w,„ and has a peak around w,„ ^ cje = lOTc as ex- 
pected. In Fig. 2(d), we show the contour plot of 
As{uJm,x). When we compare Figs. 2(a) and 2(b), we 
notice that the magnitude of the odd-frequency order 
parameter is quite small even under a strong attractive 
interaction {Vs — lOVx). Therefore, the feedback effect 
of As{uJm,x) to Ax{x) is very weak and that the value 
of Ax{x) is almost unchanged from that in the absence 
of As{ujm,x) {Vs = case). 

We show the temperature dependence of the order 
parameters in Fig. 3(a). At low temperatures, there is 
weak temperature dependence in A^, and As , while they 
are suppressed as the temperature increases. We note 
that the odd-frequency order parameter appears spon- 
taneously for r < Tc. In the vicinity of Tc, A^, show a 
square root behavior as expected. On the other hand, A^. 
show a linear temperature dependence just below Tc. To 
see this point, we also show the temperature dependence 
of A^ in Fig. 3(a). We can see that As is scaled by A^. 
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From a viewpoint of Ginzburg-Landau theory, this result 
imphes that there is a third-order term in the free energy 
for the s-wave component as 

Fs = aAi + bAi + cAsAl, (3.3) 

where a, b, and c are coefficients. We assumed real values 
for the order parameters for simplicity. In the vicinity of 
the transition temperature, both order parameters are 
small and the fourth order term can be neglected. Mini- 
mizing the free energy with respect to A^, we obtain 



2a 



(3.4) 



Since A^; shows a square root temperature dependence, 
the third-order term in the free energy explains the result 
that the s-wave component is proportional to and 
shows the linear temperature dependence. 

In Fig. 3(b), we show the coupling constant {Vs) de- 
pendence of the maximum value of |As(cjm) a^)|. We can 
see that it shows a linear dependence in the small Vg 
region. This is because the s-wave order parameter is a 
product of the induced pair amplitude and the coupling 
constant Vg. The linear Vs dependence in Fig. 3(b) im- 
plies that the constant c in Eq. (3.3) is proportional to 

Vs. 

Even though the magnitude of the odd- frequency order 
parameter is small, the odd-frequency order parameter is 
finite as shown in Fig. 3(b). The third-order term of the 
free energy given in Eq. (3.3) leads to the emergence of 
the odd-frequency order parameter below Tc of the bulk 
p-wave. 

3.2 Repulsive interaction for the odd-frequency channel 
In the previous subsections, we studied an attractive 
interaction for the odd-frequency order parameter. In the 
usual cases, superconducting order parameters are stabi- 
lized by attractive interactions. However, in the present 
case with the induced odd-frequency pair amplitude, 
the interaction is not necessary to be attractive. In this 
subsection, we consider of repulsive interaction 

(Vs < 0) and show that the odd-frequency order param- 
eter exists also in this case. This can be understood that 
the sign of the constant c in Eq. (3.4) is reversed in the 
present repulsive interaction case. To demonstrate this, 
we simply use a negative value of the coupling constant 
in Eq. (2.23) for the s-wave order parameter. 

In Fig. 4(a), we show the result for Vs = — 10T4. The 
spatial dependence of A^; (x) is basically the same as that 
for the attractive case shown in Fig. 2(a). Concerning the 
s-wave, there is no big difference between the attractive 
and repulsive cases. In Fig. 4(a), we can see that the 
sign of the order parameter is positive and is reversed 
from that of the positive interaction case [see Fig. 2(d)] 
reflecting the sign- reversed coupling constant. 

Another point is that the magnitude of the order pa- 
rameter is reduced compared to the attractive case [see 
Fig. 2(d)]. In the self-consistent calculation in the repul- 
sive case, the finite pair amplitude gives rise to a positive 
order parameter. However, the positive order parameter 
induces a negative order parameter in the next iteration 
step when we solve the gap-equation self-consistently. 
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Fig. 4. (Color online) (a) Contour plot of As{uJm,3:). The cou- 
pling constant is chosen as Vb = —lOVx. (b) Repulsive cou- 
pling constant (Vs < 0) dependence of the maximum value of 
As(wm,a;)| at a; = for T = O.lTc 



Therefore, there is a competition between the pair am- 
plitude and order parameter, where the former and lat- 
ter favor the positive and negative sign. Since the driv- 
ing force from the pair amplitude is stronger, the self- 
consistently determined order parameter becomes posi- 
tive but its magnitude is strongly reduced from that in 
the attractive case. This is also explained by the free en- 
ergy, since the a term in Eq. (3.4) becomes large in the 
repulsive interaction case. 

In Fig. 4(b), we show the repulsive coupling constant 
dependence of the maximum value of ] As(a;m, x)\ ai x — 
0. We can see that the slope becomes gradual in the 
strong repulsive coupling region which differs from the 
attractive coupling case shown in Fig. 3(b). 

This type of superconducting order parameter appear- 
ing under a repulsive interaction was previously discussed 
in a dx2^y2-wave superconductor near an interface to a 
normal metal and near a surface. ^^'^^ In the latter case, 
when a surface direction deviates from the symmetric 
axes, (1,0,0) or (1,1,0), a finite s-wave pair amplitude 
is induced owing to the symmetry mixing at the surface 
between the dx2_y2--waYe and s-wave. In this case, the in- 
duced pair amplitude has an even-frequency dependence 
and it is stabilized as an s-wave superconducting order 
parameter under a finite interaction. It does not mat- 
ter if the interaction is attractive or repulsive. ^^'^^ We 
note that the physics of the emergent order parameter 
under a repulsive interaction is essentially the same in 
the dx2^y2-wave and the present p-wave cases. 

The repulsive interaction used here is just a theoret- 
ical model. In a realistic case, a short-range Coulomb 
repulsion can be the candidate for the strong repulsive 
interaction. 
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Fig. 5. (Color online) (a) Spatial dependence of the p-wave order 
parameters at T = O.lTc for Vx = Vy. The parameter is chosen as 
u)c = lOOTc. (b) Spatial dependence of the s-wavc order param- 
eter. The Matsubara frequency is fixed as uJm = 17.3Tc at which 
|As(ti;m, has the maximum value. The set of parameters are 
chosen as = lOTc and Vs = WVx- (c) Matsubara frequency 
dependence of the odd-frequency spin triplet s-wave order pa- 
rameter at various positions, (d) Contour plot of As(u)m,,x)- 



3.3 Px + ^Py -wave case 

Next, we study the case of the chiral Px + ipy-wave 
case with Vx = Vy. The transition temperature and 
the solution of the gap equation Eq. (2.23) at T = in 
the bulk {x oo) are given by^^'^^ 

2e'''a;c 



Tr 



Ax{oo) ~ Aj,(oo) ~ 2lUcC 



1.76Tc = Ao. (3.5) 




The coherence length is defined by ^ = vf/Aq. We note 
that Aq for the px + ipy-wave is smaller than that for the 



domain wall 



Fig. 6. (Color online) Domain wall between — ipy and + 
ipy pairing states. The wavevector is conserved when the quasi- 
particle passes through the domain wall (fcpj = fcp2)- 



pure pa;-wave case studied in the previous subsection. 
Accordingly, ^ for the px + ipy-wave is larger than that 
for the px-wave. 

We next show the self-consistently determined order 
parameters in Fig. 5. The px-'wa.ve component Ax{x) de- 
creases near the surface, while the p^-wave component 
Ay{x) is enhanced as expected [see Fig. 5(a)]. As in 
the Px-wave case, the odd-frequency order parameter ap- 
pears near the surface [ see Fig. 5(b)]. It shows the same 
behavior as that in the px-wavc case as shown in Figs. 
5(c) and 5(d). 

3.4 Domain wall (px — ipy\px + ^Py) case 

In case of the Px ± ip^-wave superconductors, the 
Px ± ipy states degenerate and break the time-reversal 
symmetry. In this case, we can expect a domain wall 
formed between the two paring states. There are two 
types of the domain wall. One is formed between px — ipy 
and Px +ipy, while the other is formed between —px + ^Py 
and —px + ipy ■ We note that the former domain wall is 
energetically favorable and is stabilized against the lat- 
ter case.^^"^^ Near such domain wall, bound states exist 
owing to the sign change of the order parameter of the 
Py-wave after passing through the domain wall. In this 
case, we can also expect that the odd-frequency order 
parameter appears. 

We consider a domain wall shown in Fig. 6. We as- 
sume that there is no scattering at the domain wall and 
that the quasi-particle travels along a straight line. Since 
the odd-frequency amplitude is induced by the sign- 
changing component (py-wave), the odd- frequency order 
parameter is pure imaginary. In case of the domain wall 
Px — Wy \Px + iPy ; the induced odd- frequency amplitude 
has a spin-triplet c?a;j,-wave symmetry. We will discuss 
this point later in §4.3. Then, the order parameters are 
expressed as 

Ai(fcF,Wm,x) 

= Ax{x)cosek + iAy{x)sin9k+iAd {ojm , x) sin 29k , 
A2(fcF,a;m,x) 

= Ax{x)cosek-iAy{x)sin9k+iAd {ojm , x) sin 29k ■ 

(3.6) 

Here, iArf(w,„, a;) sin 20^ represents the d^y-wave order 
parameter. To represent the dxy-wavc symmetry, we use 
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sin20fc function for simplicity. Wc can solve the quasi- 
classical Green's function along the quasi-classical tra- 
jectory as in the scattering at the surface case. In case of 
the domain wall, the gap-equation is altered as 



Ay{x) \=T 



,(3.7) 



214 cos 9kfxi9k,0J,n,x) 

X I -i2VySin0kfy{0k,uJrn,x) 

-\2VdV- {uji,ujjn) sin26'fc/rf(6'fc, w™, x) 

where 

fx{Ok,(^,n,x) = f{dk,U!,n,x) + f* {9k,UJ„i, x) 

+ f{-6k,UJm,x) + /*(-6'fc,IU™,x), 
fy{0k,l^r,i,x) ^ f{6k,UJr,i,x) + f {0k,UJm,x) 

- f{-0k,UJm,x) - f {-0k,i^m,x), 
fd{0k,^7n,x) = f{ek,UJ,n,x) ~ f*{0k,UJ,n,x) 

- f{-0k,UJm,x) + f*{-0k,UJ,n,x).{3.8) 

The following relations were used to derive the gap equa- 
tion in the present case: 



f{0k + TT, -UJm, x) = -fi0k, iOm, x) , 
f{0k, -UJm,x) = f*{0k,UJ7n,x). 



(3.9) 



To demonstrate the emergent odd-frequency dxy-'wsive 
order parameter, we simply use V-{LUi,u!m.) defined by 
Eq. (2.26) as the frequency dependence of the effective 
interaction for the d^y-'wa.ve. 

In Fig. 7, we show the result. In the domain wall case, 
A.y{x) is suppressed near the boundary, while A.^{x) is 
enhanced as shown in Fig. 7(a). As in the surface case, 
a finite d^y-wavc order parameter appears with an odd- 
frequency dependence as shown in Figs. 7(b) and 7(c). 
In Fig. 7(d), we show the contour plot of Ad{uJm, x). 

4. Emergent Odd-Frequency Order Parameter 
and Mid- Gap Bound State 

It is well known that surface breaks the unconventional 
superconducting order parameter and induces surface 
bound state. For instance, in Copper oxide high temper- 
ature superconductors, (1,1,0) surface breaks the dx^-y2- 
wave order parameter and induces the mid-gap surface 
bound state. ^^^^^ When the time-reversal symmetry is 
broken, the energy of the mid-gap (zero-energy) bound 
state is shifted to a finite energy and it carries a surface 
current to induce a spontaneous magnetic field localized 
in the vicinity of the surfacc.^^'^^'^^'^^"^^ 

It was previously reported that the induced odd- 
frequency pair amplitude is related closely to the exis- 
tence of the bound state. ^^"^"^ In this section, we examine 
this point in terms of the quasi-particle wavefunction by 
solving the Andreev equation. 
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Fig. 7. (Color online) (a) Spatial dependence of the p-wave order 
parameters at T = O.lTc for Vx = Vy. The parameter is chosen as 
oje = lOTc. (b) Spatial dependence of the daj^-wave order param- 
eter. The Matsubara frequency is fixed as uJm = 17.3Tc at which 
\^d{^m,x)\ has the maximum value. The set of parameters are 
chosen as tjg = lOTc and Vd = lOVx- (c) Matsubara frequency 
dependence of the odd-frequency spin-triplet dxy-viaMe order pa- 
rameter at various positions. We note that /^^{uJm, is an even 
function with respect to x. (d) Contour plot of A^(a;m,a;). 



^.1 Quasi-classical wavefunction 
4-1.1 Andreev equation 

Let us begin with the following Andreev equation for 
unconventional superconductors:^^' ''^ 

/-it)F • V A{k-p,r)\ fuEirf 
l^A(fcF,r) ivp-V ) \VE{r)^ 

Here, we focus on the p^j-wave state with a real value 
of the order parameter. The Andreev equation enables 
us to investigate the quasi-particle states in terms of 
the slowly- varying function of the order of the supercon- 



= E f"^;^) ) . (4.1) 

\vE{r)> 
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ducting coherence length, where the rapid oscillation of 
the order of the Fermi wave length is averaged out. The 
derivation of the Andreev equation is given in Appendix 
A. There are positive {E > 0) and negative (E < 0) 
solutions in Eq. (4.1). They correspond to the particle 
and hole solutions, respectively, and are not independent 
each other. There is the following relation between these 
solutions: 



u-E{r)\ ^ f-VE{r) 
v-eM J \uE{r) 



(4.2) 



Summing all energy eigenvalues, the field operators for a 
fixed fcp are written as 



70 



E 

E>0 



UE 
VE 



-VE 
UE 



(4.3) 



where the definition of the field operators are given by 
Eq. (A-9). In Eq. (4.3), the first term is for the mid-gap 
(zero-energy) state. The second and third terms are for 
the particle and hole solutions, respectively. We intro- 
duced two kinds of fermions, jEf and jei- In case of the 
mid-gap state, however, we note that there is only one 
energy eigenstate for each fcp • For a fixed fcp , the Hamil- 
tonian is expressed in the following diagonal form with 
the 7 fermions for superconducting quasi-particles: 



(4.4) 



£;>o 



4-1-2 Solution of Andreev equation 

In this subsection, we focus on the Px-wave supercon- 
ductor with a specular surface at a; = 0. Since the spatial 
dependence of the order parameter is irrelevant to un- 
derstand the mid-gap state, we assume a constant order 
parameter to solve the Andreev equation. As discussed 
in Appendix B, we can solve the Andreev equation along 
the quasi-classical trajectory as shown in Fig. B-1. It is 
convenient to introduce the following functions: 



f±{r) = u{r) ± iv{r). 



(4.5) 



Here, r is the coordinate along the trajectory. We 
dropped the subscript 'i?' in u, v, and f± for convenience. 
The Andreev equation is then expressed as^^ 

Ef±{r) = -ivF^/T(r)±iA(fcF,r)/^(r). (4.6) 

We differentiate Eq. (4.6) with respect to r and obtain^^ 



E^ - A2(fcF,r) ± i;F^A(fcF, r) 
or 



(4.7) 



Along the quasi-classical trajectory, the order parame- 
ter is uniform, however, its sign suddenly changes from 
negative to positive at r = [see Fig. B-l(b)]. There- 
fore, the derivative of the order parameter is replaced as 
dA^kp, r)/dr = 2\Akp\S{r). Here, Afcp is the uniform or- 
der parameter for fcF- S{r) is the Dirac's delta- function. 
Introducing the coherence length ^ = WF/|Afcp|, we can 



rewrite Eq. (4.7) as 



f±{r)^{e'-l)f±ir). (4. 



Here, e is a dimensionless energy defined by £ = 
£'/|Afep|. We note that Eq. (4.8) has the same form of 
a Schrodinger equation with a delta-function potential. 
Since the potential is attractive, there is one bound state 
at e = for f+{r), while there is no bound state for 
/- (r) . The bound state solution is given by 



2 _M 



(4.9) 



In case of the continuum state (|e| > 1), we classify the 
solution into two types following Takayama et alJ^ 

.fi'\-r) = fi'\r), /i^VO - ^f^-\r), 

f+\-r) = -f+\r), /i'^(-r) = /l'\r). (4.10) 

The explicit form of the solution is given by 



/l^^(r)=i 



L 1 + i^fc 



2 l-iik 




{^k cos kr — sin fc|r|) . 



sin kr, 



2 /1-i^fc . , 

sm kr, 



LV 1 + i^k 



Ll + i^k 



(^fc cos fcr -I- sin fc|r|) , (4.11) 



where the energy eigenvalue is given by e = -^/l -f (C^)^. 
The positive wave numbers fc > are sufficient for 
the solution, since the wavefunctions are identical under 
k — >■ — fc except for the phase factor. The normalization 
condition for f± is chosen as''^'''^ 



[|/+(r)|2 + |/_(r)|2] =^ 

The corresponding u and v are given by 

Voir) 



(4.12) 



, , 1 _^ 



2^ 



1 _M 

i^=e 5 , 



u(i)(r) = -xl- 



1 



2 V L l-f i^fc 



^k cos kr' — sin fcjr'l + i\/ 1 + (<^fc)^ sin kr' 
I 



.1 2 



2 V i 1 + i^fc 



X ^fc cos kr' — sin fcjr'l — i-\/ 1 + (^fc)^ sin kr' 

(2)r N 1 /2" 1 
'(r) = -\ 



i^k 



^fccosfcr' + sinfc|r'| + (^fc)2sinfcr' 



1 



2 \ Ll + i£_k 
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COS kr' + sin k\r'\ — i-\/l + (^/c)^ sin kr' . 

(4.13) 

Here, uq and vq are wavefunctions for the mid-gap bound 
state, while u^*^ and (i = 1, 2) are for the continuum 
states. We note that u and v satisfy the following nor- 
malization condition: 



dr [|u(r)p + |i;(r)p] = 1. 



(4.14) 



4.2 Superconducting pair amplitude 
4.2.1 General description 

The superconducting pair amplitude is defined by the 
following form: 

= -{,p^{r,T)i;^{r'))0{T) + {^P^{r')^b^{r,T))e{-T). 

(4.15) 

Here, Tt represents the time-order operator with respect 
to the imaginary time t and 9{t) is the step function. 
Substituting Eq. (A-3) into Eq. (4.15), we obtain 

Ft4.(r,r',r) 

= _ ^ e'^--e-"=-'-'(T,^fc,t(r-,r)V-fcp'i(r')) 



feF 



X 1- 



2 

X 



(4.16) 



where i? and x represent the center of mass and relative 
coordinate of the Cooper pair, respectively. Fkp-fi{r,T) 
is defined by 

Fk^nir, r) = - {Tr^PkMr, r)V-fcF;(r-)) (4.17) 

as the Fourier transformed pair amplitude. We note that 
r in Fkp^i{r,T) corresponds to the center of mass coor- 
dinate of the Cooper pair. The Fourier transformed pair 
amplitude with respect to r is defined by 



■FfcFti('''i^™) 



dTe'""^^Fk,nir,r). (4.18) 



Here, uJm is the fcrmionic Matsubara frequency and /? = 
1/T. 

In the following discussion, we express Fkj^^i{r,iuJm) 
with the energy eigenstates of the Andreev equation and 
examine what state forms the odd-frequency pair ampli- 
tude. For this purpose, we express the pair amplitude by 
using the wavefunctions u and v. Substituting Eq. (4.3) 
into Eq. (4.18), we obtain 



r/3 



dre" 



/ 

^0 



dre" 



{uo(r)wo(r)(7o7,5) 



E>0 



r^«o(r)t;S(r) + ^2Re 



E 



(4.19) 



Here we assumed a low temperature limit. The first term 
is from the mid-gap bound state, while the second term 
is from the continuum states. It is clear that the first 
term has an odd-frequency dependence. Equation (4.19) 
clearly indicates that the mid-gap bound state generates 
the odd-frequency pair amplitude. It is the key to under- 
stand the emergent odd-frequency order parameter. 

4.2.2 Pair amplitude near the boundary 

We discuss next the pair amplitude in terms of the so- 
lution of the Andreev equation. Substituting Eq. (4.13) 
into Eq. (4.19) and summing up the two types of solu- 
tions, we can divide the pair amplitude into the even- 
and odd-frequency components as 

Fk.ni^, i^™) - F^lf^ir, iw„0 + F^^^^ir, iu;^). (4.20) 

Here, F^'^°^{r,iuJm) and F^^^j^(r, iw^) are the even- and 
odd-frequency components defined by 



TTicvcn 



(r, iuJn 



|Afep|sgn(r) 



k>0 



F^:U^.icOra) 



1 



lAfcF 



-e i 



+ + [vFky 



vpk 



■ sin kr, 



fc>0 



■ sin 2fc|r|. 
(4.21) 



The even-frequency component is an odd function with 
respect to r refiecting the symmetry of the p^^-wave. We 
note that only the continuum states contribute to the 
even- frequency component. On the other hand, the odd- 
frequency component is an even function with respect 
to r. This indicates that it has an s-wave symmetry. In 
the odd-frequency component, both of the bound and 
the continuum states contribute to the pair amplitude. 
In the following subsections, we study the two frequency 
components separately. 

4.2.3 Odd- frequency pair amplitude 

We first replace the summation over k with its integral. 
The pair amplitude is then expressed as 



1 

TT 



dfc- 



v-pk 







UJr, 



sin2fc|r| 



. 1 r , 1 

I7; — / dz: 

2WF J -00 ^^rn - Z 

X ^— Im 

— TT 
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. 1 |AfcJ 2,/^„^+A|Jr|/t,F 



(4.22) 



+ ivpky and used the 



where we introduced z - 
Lehmann representation with an infinitesimal small pos- 
itive number S. In terms of the x{> 0) coordinate, where 
X = \vpx\\f\/vF, we obtain 



(4.23) 



We can see that the odd-frequency pair amplitude has 
the s-wave symmetry, i.e., iw™) is symmetric 

with respect to /cfx ^ — ^Fa; and fcpy ^ —kpy The 
pair amplitude decreases exponentially as x increases, 
reflecting the fact that it is formed mainly by the mid-gap 
bound state. We note that Eq. (4.23) is consistent with 
the result of the odd-frequency pair amplitude obtained 
by the quasi-classical Green's function {gi component) 
given in Eq. (2.2). 

4.2.4 Even-frequency pair amplitude 

The even-frequency pair amplitude is rewritten as 

1 



dk- 



_ 1 
X ^— Im 

— TT 



TT JO 
1 



^#^^T^2sin^fc. 

+- A2 + («Ffc)2 



|Afcp|sgn(r) 



^ 1 -|AfcJsgn(r) / _ ^-2^a;„2 + A^^|r|/. 



2^f,/W7Ta 



(4.24) 



In terms of the x coordinate, where |Afcp|sgn(r) — Aj-p 
we obtain 



1 -Afe, 



2VF 



A? 



1-e 



(4.25) 



We can see that the even-frequency pair amplitude has 
the px-wave symmetry. We note that Eq. (4.25) is consis- 
tent with the result of the even- frequency pair amplitude 
obtained by the quasi-classical Green's function (52 com- 
ponent) given in Eq. (2.2). 

4.2.5 px + ipy-wave case 

In the chiral Px + ip^-wave case, the bound state en- 
ergy is located at = Aj,(fcF) for each fepy^^'®* and 
there is no mid-gap bound state except for Ay{kp) = 0. 
For kpy > 0, the bound state energy is positive, while 
it is negative for kpy < 0. Thus, the bound states are 
not symmetric with respect to kpy and they carry a fi- 
nite surface current reflecting the broken time-reversal 



symmetry.^^'^^ 

When the time-reversal symmetry is broken, the rela- 
tion given in Eq. (4.2) does not hold and we have to solve 
both the particle and hole solutions. Let us consider the 
bound state contribution. In the present px + ipy-wave 
case, the bound state contribution in Eq. (4.19) is re- 
placed as 



" W]*, (4.26) 



-■ Uo[r)Voir) ^ -. r— 

lUJ,n lUJrn-Ay{kF) 

where u^{r) and v^{r) are wavefunctions for the bound 
state at E = Aj^(A;f). For constant Axikp) and A.y(fcF), 
they are given by 



lAxikp)] i^x(fcF)i 



2'yF 



|A,(fcF)| 



Ax(fcF)l 



|r| 



2VF 



(4.27) 



Substituting Eq. (4.27) into Eq. (4.26), we obtain 
1 



iu},n - Ayikp) 

1 |A,(fcF 



u^ir)[v^ir)Y 



A2(fcF 



iA.y(fcF)]e 



(4.28) 



Here, the Um term in the numerator represents the odd- 
frequency amplitude, while the — iAy(fcF) term is for the 
even-frequency amplitude. When we add the contribu- 
tion from the continuum states, the pair amplitudes are 
expressed as^^ 



F^^'^^{x,iu^) 
FZ'^^l{x^Lu^) 



1 bJm\Ax{kY)\ ^_ 

2uF a;™2_|. A2(feF) 
1 -A,(fcF 



2qx 



2UF ^ 

\Ay{kY) 



where 



2VY 



1 



(l-e 
A^(fcF) 



2qx\ 



0J„ 



A?(fcF 



A2(fcF) + A2(fcF), 



fl/\vFx 



(4.29) 



(4.30) 



It is clear that the odd-frequency pair amplitude exists 
also in the chiral px + ipy-wave state. As in the px- 
wave case, the odd-frequency pair amplitude decreases 
exponentially with x. The first and second terms in 
-^fept" '^'^■m) represent the Px- and py-wave components, 
respectively. We can see that the former vanishes at 
X = 0, while the latter is enhanced. It is consistent with 
the spatial dependence of the order parameters shown in 
Fig. 5(a). 

4.3 Symmetry of odd-frequency order parameter 

In this subsection, we discuss the symmetry of the 
emergent odd-frequency order parameter. As in Eqs. 
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71 



(b) 



+ 



Fig. 8. (Color online) Order parameter [thin (red) lines] along 
the quasi-classical trajectory. The solid (blue) lines represent the 
odd-frequency order parameter near the boundary. The ± sign 
represents the sign of the order parameter, (a) The case where 
the sign of the order parameter changes from negative to positive, 
(b) From positive to negative case. 



(4.19) and (4.26), the odd- frequency pair amplitude is 
generated by the bound state. In the presence of a fi- 
nite interaction for the odd-frequency component, the 
pair amphtude is stabilized as a superconducting or- 
der parameter. Since the bound state is formed by pair- 
breaking perturbations such as surface and domain wall, 
we can understand that the odd-frequency order param- 
eter is brought about by the pair-breaking effects. In case 
of the surface and domain wall, the spin-triplet nature of 
the p-wave Cooper pair holds when the quasi-particle is 
scattered by the surface or when it passes through the 
domain wall. This means that the induced order param- 
eter has the same spin-triplet property. Concerning the 
symmetry of the orbital of the Cooper pair, we saw in 
the previous section that the s-wave and (ia;j,-wave order 
parameters were stabilized near the surface and domain 
wall, respectively. This means that the orbital symmetry 
of the odd-frequency order parameter depends on the ge- 
ometry of the boundary relative to the symmetry of the 
unconventional order parameter. We discuss this point 
below. 

We first summarize the result in the previous subsec- 
tions. When the sign of the order parameter changes from 
negative to positive along the quasi-classical trajectory, 
a positive pair amplitude is induced. For an attractive 
interaction, the sign of the odd-frequency order param- 
eter becomes negative as shown in Fig. 8(a). In the op- 
posite case, where the order parameter changes its sign 
from positive to negative, the sign of the order param- 
eter becomes positive [see Fig. 8(b)]. This is the key to 
understand the orbital symmetry of the odd-frequency 
order parameter. 

Let us begin with the p^-'wave superconductor. At the 
surface, there are two types of scatterings. One is shown 
in Fig. 9(a), while in the other one the directions of 
the wave vectors are reversed. In both cases, the sign 
of the order parameter changes from negative to positive 
as shown in Fig. 8(a). Therefore, the odd-frequency or- 
der parameter has an s-wave symmetry with a negative 
sign as shown in Fig. 9(b). This is consistent with the 
results in Fig. 2. This property holds also in the Px + '^Py- 
wave case, since only the p^^-wave component carries the 
sign-change of the order parameter. 

In case of the px — ^Py\Px + ^Py domain wall, the quasi- 
particle passes through the domain wall. The sign-change 
of the order parameter occurs in the Py-wave component 
in this case. For fcpi shown in Fig. 9(c), the sign changes 
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Fig. 9. (Color online) (a) Incident and reflected quasi-particles at 
the surface, (b) Orbital symmetry (s-wave) of the odd-frequency 
order parameter, (c) Incident and reflected quasi-particlcs for 
the domain wall, (d) Orbital symmetry (d2;y-wave) of the odd- 
frequency order parameter. 



from negative to positive, while it is opposite for h^i. 
The former and the latter cases correspond to the sit- 
uation shown in Figs. 8(a) and 8(b), respectively. This 
leads to a d^y-wavc symmetry of the odd-frequency order 
parameter with a negative sign of the global phase factor 
(tt phase shift relative to the sign-change even-frequency 
component) as shown in Fig. 9(d). This is also consistent 
with the results shown in Fig. 7. Since the sign-change py- 
wave component is pure imaginary for \hsi px~\-Py\px^'vpy 
domain wall, we note that the odd-frequency dxy-wakve 
order parameter becomes also pure imaginary. 

5. Summary and Discussions 

In this paper, we investigated the emergent odd- 
frequency superconducting order parameter near bound- 
aries in unconventional superconductors. We focus on 
the two-dimensional spin-triplet even-frequency Px- and 
Px + ipy-wave superconducting states of the d \\ z type, 
where the latter is the candidate of the paring state for 
Sr2Ru04 superconductor. 

In the px-wave superconductors with a specular sur- 
face along the y-direction, it is known that the order 
parameter is suppressed and the mid-gap bound state 
appears near the surface. In the previous studies, it was 
suggested that the odd-frequency spin-triplet s-wave pair 
amplitude is induced near the surface. In our study, 
we took the idea one step further and investigated the 
odd-frequency order parameter by introducing a finite 
interaction in the odd-frequency pairing channel. In par- 
ticular, we considered the electron-phonon interaction, 
since it leads to a finite attractive interaction in the spin- 
triplet s-wave channel. We extended the quasi-classical 
theory to include the frequency-dependence of the order 
parameter and confirmed that the conventional relation 
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between the particle and hole components of the order 
parameter, i.e. Eq. (2.10), must be used when we con- 
sider the induced odd-frequency order parameter with 
the primary even-frequency order parameter. Solving the 
extended Eilenberger equation self-consistently, we de- 
termined both the frequency and spatial dependences of 
the order parameter and found that the odd-frequency 
spin-triplet s-wave superconducting order parameter ap- 
pears near the surface (see Figs. 2 and 5). 

To understand the emergence of the odd-frequency or- 
der parameter, we expressed the pair amplitude with the 
quasi-particle wavefunctions as given by Eq. (4.19). It 
clearly indicates that the mid-gap bound state generates 
the odd-frequency pair amplitude. In case of the chiral 
Px + ipy-wavc superconductor, the bound state energy is 
located at E = Ay{kp) for each fcpj, reflecting the bro- 
ken time-reversal symmetry. Although they are not the 
mid-gap state, the surface bound state generates the odd- 
frequency pair amplitude. Under a finite interaction, the 
pair amplitude is stabilized and the odd-frequency order 
parameter appears. 

When we compare the magnitude of the order param- 
eters, we notice that the odd-frequency component is 
quite small even under a strong attractive interaction 
in the odd-frequency channel. The main reason for this 
is that the odd-frequency order parameter is hard to be 
stabilized in the bulk owing to the limiting frequency 
region for the attractive interaction.^^ This property re- 
mains also in the non-uniform case of the induced odd- 
frequency order parameter near boundaries. The other 
reason is that the frequency dependence in the bulk p- 
wave order parameter was not taken into account un- 
der assuming the weak coupling theory. In this case, 
the p-wave component is overestimated and dominates 
the minority s-wave odd-frequency component. Although 
the magnitude of the odd-frequency order parameter is 
small, we elucidated how it appears near boundaries in 
unconventional superconductors. It appears as long as 
the coupling is finite irrespective of its sign (attractive 
or repulsive) as shown in Figs. 3 and 4. This point can 
be explained qualitatively by the third-order term in the 
free energy given in Eq. (3.3). When the odd- frequency 
order parameter is finite, physical quantities show differ- 
ent behavior. Although the effect is expected to be weak 
in case of a small odd-frequency order parameter, it is 
important to study this point in a general case. 

In this paper, we focussed on surface and domain wall 
in unconventional superconductors. The idea of the emer- 
gent odd-frequency order parameter can be adopted also 
to interfaces between normal metal and unconventional 
superconductors or between ferromagnet and supercon- 
ductors and so on, since the odd-frequency pair ampli- 
tude is induced in these cases. ^^'^"^ In particular, in the 
former larger magnitude of the odd-frequency 

order parameter is expected in the normal metal side, 
since the unconventional (even-frequency) order param- 
eter vanishes in the normal metal and there is no compe- 
tition between the even- and odd-frequency ones. In this 
case, we can expect that the effect of the odd-frequency 
order parameter is seen in the normal metal side. These 
points will be examined in the future work. 
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Appendix A: Derivation of Andreev Equation 

In case of the px-wa.ve of the d-vector parallel to the 
z-axis type, the field operators satisfy the following equa- 
tion of motion: 



dr V^j(r,r) 



dr' 



(5(r - r')£(-iV') 
-A*(r,r') 



A(r,r') 



-6{r - r')e{-\\/'] 



(A-1) 

Here, ip^{r,T) and 'i/'j (r,r) are field operators in the 
Heisenberg representation defined by 



V'l(r,T) 



V^l(r) 



-Ht 



(A-2) 



where T-L is the Hamiltonian. A(r,r') is the position de- 
pendent order parameter. e(— iV') = (— iV')^/(2m) — i^p 
is the operator for the kinetic energy measured from the 
Fermi energy E-p. 

To derive the Andreev equation, we rewrite the field 
operators as 



V't(r) 



(A-3) 



The field operators ip^{r) and ip'^ir) contain the rapid 
oscillation term e''"'' '', while the field operators tpk-pti''') 
and i^-kpii''') describe slowly varying component of the 
order of the superconducting coherence length. We sub- 
stitute Eq. (A-3) into Eq. (A-1) and apply the quasi- 
classical approximation,'*^'™ where the relative coordi- 
nate X = r r' and the center of mass coordinate 
R ~ {r + r')/2 arc introduced. The two coordinates, 
r and r', are then expressed asr ~ R + x/2 and 
r' = R x/2. The order parameter is written by x and 
R as 



A(r,r') A{x,R) 



(A-4) 



The order parameter and the slowly varying field opera- 
tor can be expanded as^^ 

A(a;, R) ^ A(a;, r - ^) ~ A(a;, i')-^- VA(a;, r) H , 

^kA'^') = ^kA'^ - a;) ~ ^kAf) - X ■ VV^fcp(r) H . 

(A-5) 

For a fixed fcp, we substitute Eq. (A-5) into Eq. (A-1) 
and obtain^* 

J dr'A{r,r')i>lir') ^ J dr'A(a;, fi)V'Ip^(r')e"^'' ''' 
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e 



dx 



A{x,r) - - • VA(a;,r) 



= e"=-'^A(fcF,r)^fc,(r), 
where we introduced 



(A-6) 



A(fcF, r) = J dxA{x, r)c-''^''-^ (A-7) 

as the Fourier transformed order parameter with respect 
to the relative coordinate of the Cooper pair. In A(fcF, r), 
the fcp dependence represents the orbital symmetry of 
the Cooper pair, while the slowly varying r dependence 
represents the spatial dependence of the order parame- 
ter in terms of the center of mass coordinate. Then, we 
obtain the following equation of motion of the field op- 
erators:^^ 

-it>F-V A{kF,r)\ f ipkp^{r,T) 
A{kp,r) ivF-W J [^l^^^ir,T] 



9t W-k^i{r,T) 



,(A.8) 



where we used — A*(— fcp, r) — A{kp,r) for the p^j-wave 
with a real value, = fep /m represents the Fermi veloc- 
ity. To find energy eigenstates, we put the field operators 
as 



i^k^^ir) \ _ fuE{r)\ 



(A-9) 



Here, is a fermion operator for an energy eigenvalue 
E. We assume that the Hamiltonian has a diagonalized 
form as E'^I^'^e- Substituting Eq. (A-9) into Eq. (A-8), 
we obtain the following Andreev equation: 

-ivF ■ V A{kF,r)\ / U£;(r)\ ^ ^ f UE{r) 
A{kF,r) it^F-V J \VEir) J " \VEir) 

(A-IO) 

Appendix B: Boundary Condition of Andreev 
Equation 

In this appendix, we discuss the boundary condition of 
the Andreev equation at the surface. The quasi-classical 
trajectory is shown in Fig. B-l(a). We take the coordinate 
r along the quasi-classical trajectory. The incident quasi- 
particle with fcpi wavevector comes from r = —oo and 
is reflected at the surface at r = 0. The reflected quasi- 
particle travels with k-p2 wavevector toward r = oo. The 
order parameter for the incident quasi-particle is A^p^, 
while it is A^p^ for the reflected quasi-particle. For an 
uniform order parameter, the spatial dependence of the 
order parameter along the quasi-classical trajectory sud- 
denly changes its sign at the surface as shown in Fig. 
B-l(b). We solve the Andreev equation (A-10) in the two 
regions separately and obtain the solutions analytically. 



r= 




r> 






X 













r< 







(b) 



Fig. B-1. (Color online) (a) Quasi-classical trajectory, (b) Spa- 
tial dependence of the order parameter along the quasi-classical 
trajectory. 



Let us discuss the boundary condition of the wavefunc- 
tion of the quasi-particle. In the presence of the specular 
surface along the j/-direction, the y component of the 
wavevector {k-py) is conserved. The y dependence of the 
wavefunction is then expressed by c^^^^yy . On the other 
hand, the ifep^; components are coupled by the surface. 
For a flxed fcpy, the field operator is written as 



(B-1) 



Here, fcpi = {-kYx^kpy) and k-p2 = {kYx,kYy). The y 
dependence is extracted as e''^'^"^. The field operators 
V'fcFt(^) '/'-fcF^(^) ^"^^-^ have the x dependence. Since 
the field operator must vanish at the surface, we obtain 



V'fcFiTW + ^fcF.rW-o, 
V'U.;(o) + '^-fcF.;(o) = 0' 



(B-2) 



for each k-py. These lead to the following boundary con- 
dition for u and u at r = 0: 



u(0)^ 



-(0)/ .F 



(B-3) 



Since the Andreev equation (A-10) is invariant under the 
change of the global phase of the wavefunctions, the mi- 
nus sign in Eq. (B-3) can be absorbed in the wavefunc- 
tions. Therefore, we solve the Andreev equation in both 
r < and r > regions separately and connect the so- 
lutions with the following conditions at r ~ 0: 



"(0) 
v{Q) 



v{0) 



(B-4) 



This indicates that it is sufficient to solve the Andreev 
equation (A-10) continuously along the quasi-classical 
trajectory. In the similar way, the quasi-classical Green's 
function is also solved continuously along the quasi- 
classical trajectory as we mentioned in Eq. (2.18). 
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